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We analyze the quantum phase transition between a semimetal and a superfluid in a model of attractively 
interacting fermions with a linear dispersion. The quantum critical properties of this model cannot be treated 
by the Hertz-Millis approach since integrating out the fermions leads to a singular Landau-Ginzburg order 
parameter functional. We therefore derive and solve coupled renormalization group equations for the fermionic 
degrees of freedom and the bosonic order parameter fluctuations. In two spatial dimensions, fermions and 
bosons acquire anomalous scaling dimensions at the quantum critical point, associated with non-Fermi liquid 
behavior and non-Gaussian order parameter fluctuations. 
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I. INTRODUCTION 



As noticed already more than 30 years ago by Hertz, ^ quan- 
tum phase transitions in correlated fermion systems may be 
described in terms of the order parameter field alone when one 
integrates out fermions from the path integral in one stroke 
and subsequently deals with an exclusively bosonic theory. 
The resultant eff'ective action is expanded in powers of the 
bosonic field and often truncated after the 0^-term. Al- 
though an analysis in terms of 0^-type theories seems mun- 
dane, the presence of two relevant energy scales, one given by 
temperature and the other given by a non-thermal control pa- 
rameter acting as a mass term in the bosonic propagator, gives 
rise to a rich finite temperature phase diagram.^ 

The Hertz-Millis approach relies on integrating out 
fermions first. In general, however, the fermions are gapless 
and consequently may lead to singular coefficients in the ef- 
fective bosonic action.«^ In such cases, it is advisable to keep 
the fermions in the theory and treat them on equal footing 
with the bosons. Various coupled Fermi-Bose systems ex- 
hibiting quantum criticality have been analyzed previously by 
(resummed) perturbation theory. ^'^'^'^'^ The complicated inter- 
play of two singular propagators promotes a controlled pertur- 
bative treatment to a formidable task. 

A suitable tool to cope with the interplay of fermionic and 
bosonic fluctuations is the functional renormalization group 
(RG) formulated for fermionic and bosonic fields. ^^'^^'^^'^^ 
The mutual feedback of gapless fermions coupled to mass- 
less bosons has been studied already with functional flow 
equations for quantum electrodynamics,^^ non-abelian gauge 
theories, and the Gross-Neveu model. 

The purpose of the present paper is to show how one can - 
within the coupled Fermi-Bose RG framework - obtain non- 
Fermi liquid properties of quantum critical fermion systems 
for which the Hertz-Millis approach is not applicable. To 
this end, we study a "Dirac cone" model of attractively in- 
teracting fermions with a linear dispersion, which exhibits a 
quantum phase transition from a semimetal to a superfluid at 
a finite interaction strength as shown in Fig. [T] The coupled 
fermion-boson flow equations are physically transparent and 
a relatively simple truncation seems to capture the essential 
renormalizations. As a central result, the fermion and order 




FIG. 1: (Color online) Schematic phase diagram of the attractive 
Dirac cone model in the plane spanned by the inverse interaction and 
temperature. The quantum critical point (QCP) at a critical interac- 
tion strength Uc separates the semimetal from the superfluid. In this 
paper, we approach the QCP at T = from the semimetallic phase 
as indicated by the arrow. 



parameter propagators develop anomalous power law depen- 
dences on frequency and momentum at the quantum critical 
point (QCP) in dimensions J < 3. 

In Sec. II, we introduce the attractive Dirac cone model and 
the associated Fermi-Bose action. Subsequently, in Sec. Ill, 
we show that a mean-field treatment of this model leads to a 
semimetal-to- superfluid quantum phase transition at a critical 
interaction strength Uc- In Sec. IV, the RG method, truncation 
and flow equations are presented. Results for the quantum 
critical behavior follow in Sec. V. We finally summarize and 
conclude in Sec. VI. 



II. DIRAC CONE MODEL 

We consider an attractively interacting fermion system with 
the bare action 



M= - r > 
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The variables k = (ko,k) and q = (qo,q) collect Matsub- 
ara energies and momenta, and we use the short-hand nota- 
tion = ^ J for momentum and energy integrals; 

includes also a spin sum. Our analysis is restricted to 
zero temperature such that the Matsubara energies are con- 
tinuous variables. The dispersion of the fermions is given by 
the "Dirac cone" 



= ±v/|k| , 



(2) 



corresponding to relativistic massless particles with positive 
(plus sign) and negative energy (minus sign). The chemical 
potential is chosen as yu = 0, such that in the absence of in- 
teractions states with negative energy are filled, while states 
with positive energy are empty. The Fermi surface thus con- 
sists only of one point, the "Dirac point" at k = 0, where 
the two cones in Fig. [2l intersect. The momentum integra- 
tion in Eq. ([T]) and all other fermionic momentum integrals 
in this article include the sum over positive and negative en- 
ergy branches of ^k- Momentum integrations are cut off' in the 
ultraviolet by the condition |^kl < Aq. 

The dispersion Eq. ([S]) is a simplified version of the dis- 
persion for electrons moving on a honeycomb lattice as in 
graphene, where the momentum dependence is entangled with 
a pseudospin degree of freedom related to the two-atom struc- 
ture of the unit cell.^^ 

For attractive interactions, the coupling constant U is neg- 
ative and drives spin singlet pairing associated with sponta- 
neous breaking of the global U(l) gauge symmetry. There- 
fore, we decouple the Hubbard interaction in the s-wave spin- 
singlet pairing channel by introducing a complex bosonic 
Hubbard-Stratonovich field cp conjugate to the bilinear com- 
posite of fermionic fields U fj^i//k+i^i//_k+ii}^ This yields a 
functional integral over i//, ij/ and with the new bare action 

Toil//, {lr,(p]= - f ifrkAiko - ^k) ^/^kcr - f fqJJ^Pq 
Jkcr Jq ^ 
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where 0* is the complex conjugate of 0, while and are 
algebraically independent Grassmann variables. The boson 
mass S = l/U plays the role of the control parameter for the 
quantum phase transition. 



III. MEAN-FIELD THEORY 

Neglecting bosonic ffuctuations by replacing (pq with its ex- 
pectation value (pq^o, the saddle-point approximation solves 
the functional integral of Eq. ([3]) and leads to the standard BCS 
gap equation 
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FIG. 2: (Color online) "Dirac cone" dispersion in two dimensions. 



The critical interaction strength Uc, beyond which the gap 
equation has a solution with non-zero 0o, is given by 

-Ao 

-Ao ^0 



where the density of states has the form 
N(0 



(6) 



with Kd being defined by f -0 = Kd f d\k\ \k\^-\ Since the 
density of states vanishes at the Fermi level (for d > 1), the 
system is stable against pairing for a weak attraction. Instead, 
a finite attraction beyond a certain threshold is necessary to 
cause superffuidity even at zero temperature. Performing the 
integrations in Eq. ([5]), we obtain the mean-field position of 
the quantum phase transition in Fig. [T] In two dimensions, we 
have: 
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Upon setting v/ = Aq = 1, the numerical value for the mean- 
field control parameter value is = 0.159. 

A similar quantum phase transition from a semimetal to a 
superffuid has been discussed previously for attractively in- 
teracting electrons on a honeycomb lattice, in the context of 
ultracold atoms and graphene. 

In the following, we take ffuctuations into account by per- 
forming a renormalization group analysis which enables us to 
compute the renormalized position of the QCP as well as the 
quantum critical exponents at and in the vicinity of the QCP. 



IV. RG METHOD 

We derive ffow equations for the scale-dependent eff'ec- 
tive action F^ [0,0,0] within the functional RG framework 
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for fermionic and bosonic degrees of freedom in its one- 
particle irreducible representation. ^^'^^'^^'^^ Starting from the 
bare fermion-boson action F^^^^ [j/^, 0^, 0] = ro[</^, 0] in 
Eq. ([3]), fermionic and bosonic fluctuations are integrated 
simultaneously, proceeding from higher to lower scales as 
parametrized by the continuous flow parameter A. In the in- 
frared limit A ^ 0, the fully renormalized efl'ective action 
r^^^[0, 0, 0] is obtained. 

The flow of is governed by the exact functional flow 
equation " " ^ 



igiijZB 



-^r^[(A,0^,0] = Str 



r(2)A[0, 0, 0] -hT?^ 



(8) 



where Y^^^^ denotes the second functional derivative with re- 
spect to the fields and is the infrared regulator (to be spec- 
ified below). The supertrace (Str) traces over all indices, with 
an additional minus sign for fermionic contractions. 

When evolving F^, infinitely many terms involving 
fermionic and/or bosonic fields with possibly complicated de- 
pendences on frequencies and momenta are generated, neces- 
sitating a truncation of the eff'ective action. 



particle-particle bubble. For small momenta and frequencies 
this dependence is quadratic, leading to the following ansatz 
for the bosonic quadratic part of the action: 

Y^.^= r 0; [z,(^g^q2) + ^]0,. (11) 

Jq 

Note that there is no linear term in frequency here, as we con- 
sider the half-filled band and the usual (imaginary) linear fre- 
quency part of the particle-particle bubble vanishes exactly 
due to particle-hole symmetry. The parameter 6 controls the 
distance to the quantum phase transition and is also renormal- 
ized by fluctuations. If the initial value of d is tuned so that 
^ ^ for vanishing cutofl' A ^ 0, we are in the quantum crit- 
ical state. The boson propagator, parametrized by two flowing 
parameters, and reads 



-1 



(12) 



The initial condition for is = 0, and for the control 
parameter^ = l/U. 



A. Truncation 

We now explain how the efl'ective action is truncated with 
the objective to capture the essential renormalization efl'ects. 

1. Fermion propagator 

To account for the renormalization of the fermionic 
single-particle properties by order parameter fluctuations, the 
quadratic fermionic term in the action is multiplied by a field 
renormalization factor Zj, 

r^(A = ~ I ^kcrZfiiko - ^k) l/^kcT , (9) 

Jk(T 

yielding the fermion propagator 

Z"^ 

Gf(k) = -{i/fM = • ^1^) 

A diverging Z/ suppresses the quasi-particle weight to zero, 
leading to non-Fermi liquid behavior. The Fermi velocity is 
not renormalized separately but kept fixed, since for the linear 
dispersion relation, Eq. (|2]), the renormalization of the term 
proportional to ^0 and the one proportional to ^k in F^^ are 
essentially the same. The initial condition for Zf is Zf = 1. 

2. Boson propagator 

The bosonic quadratic part of the bare action, Eq. ([3]), 
consists only of a local mass term. Integrating out fluctua- 
tions, a momentum and frequency dependence is generated, in 
particular by fermionic contributions involving the fermionic 



3. Order parameter self interaction 

In the Hertz-Millis theory of quantum criticality^'^ one ex- 
pands the order parameter self interaction in (even) powers of 
the bosonic fields: 

Fi[0] = ^^2« r (0*0)" , (13) 

n>2 ^ 

where the coefficients U2n are generated by fermion loops with 
2n fermion propagators and 2n external bosonic legs as shown 
in Fig. [3l At finite temperatures or when the fermions are 
gapped, the integrals corresponding to these loops are finite. 
With the bare, gapless fermion propagator at zero temperature, 
however, the loops diverge for vanishing external momenta: 

where A is an infrared cutoff' introduced here merely to 
specify the degree of the divergence. In particular, the 0"^- 
interaction diverges as W4 ~ 1 / A in two dimensions. For the 
present case of fermionic loops with pairing vertices, there 
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FIG. 3: Fermion loop which generates the bosonic 2n-point vertex, 
here fom- 5. 



are no cancellations which reduce the divergence below the 
level of simple power counting. This is in stark contrast 
to the situation for instabilities driven by forward scattering, 
where oscillating integrands and cancellations between per- 
mutations of external legs render the loops finite in the low 
energy limit.i^^ 

The divergence of the fermion loops makes the Hertz-Millis 
approach of expanding the eff'ective action in powers of the or- 
dering field alone inapplicable. Besides issues related to the 
increasing degree of divergence of higher order interactions, 
already at the quartic level one would face the problem to 
parametrize the complicated singular momentum dependence 
of the interaction generated by the fermionic loop. Integrat- 
ing the fermions along with the bosons in a coupled RG flow, 
the vertices remain finite and can be parametrized by a mo- 
mentum independent running coupling as usual. In particular, 
the 0^-term, which is generated by fermionic fluctuations, but 
then also influenced by bosonic fluctuations, can be written as 

^\<P\' = ^ I fq+pfq'-p(Pq'(Pq ' (1^) 
O Jq,q',p 

The initial condition for w is w = 0. 
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FIG. 4: Feynman diagrams representing the flow equations. 



4. Fermion-boson vertex 



The fermion-boson vertex has the form 



= g \ {if^-k+'ji^k^^(Pq-^^k+'j^^-k^ifq)- (16) 
Jk,q 

The fermion-boson vertex is not renormalized within our trun- 
cation. The usual one-loop vertex correction, which is for- 
mally of order , vanishes in the normal phase due to particle 
conservation.^^ Hence, the coupling g remains invariant at its 
bare value ^ = 1 in the course of the flow. 



B. Flow equations 

In this subsection, we derive the flow equations for our trun- 
cated efl'ective action. Both propagators, Eqs. ([TQl[T2l) , display 
singularities for vanishing momenta and frequencies, the zero- 
temperature fermion propagator everywhere in the semimetal- 
lic phase, and the boson propagator at the QCP when the 
bosonic mass vanishes. These infrared singularities are regu- 
larized by adding optimized momentum cutoflfs^^ for fermions 
(subscript f) and bosons (subscript b). 



4(k) 



l^kl] 



v^v-; - Z/(-Asgn[^k]+^k)^[A- 
<(k) = Z,(-A2 + q2)^[A2-q2] , (17) 

to the inverse of the propagators. We denote the regularized 
propagators by G/r and GhR. Note that we have set v/ = 1, 
such that A is a common momentum cutofl' for fermions and 
bosons, which is a frequent choice for critical fermion-boson 
theories The scale-derivatives of the cutofl's read. 



OaR'^ =R'> = -Zf sgn [^k] ^ [A - l^kl] 



V 



R^ = -2ZtAe[A^-q^] 



where terms proportional to Zf and are neglected. These 
additional terms are of higher order in the vertices. Arguments 
buttres sing this commonly employed approximation are given 
in Ref.M 

The recipe to obtain the flow equations is now the follow- 
ing: one executes a cutofl'-derivative acting on R^^ in the an- 
alytic expressions corresponding to all 1-loop one-particle- 
irreducible Feynman diagrams for the parameters Z/, Z^, S, 
u, which can be constructed with the interaction vertices and 
the regularized propagators GfR and G^r, as shown in Fig. |4l 
For the cutofl' derivative and loop integration we use the short- 
hand notation 

The feedback of bosonic fluctuations on the fermionic prop- 
agator is captured by the self-energy diagram in the first line of 
Fig.m leading to the flow equation for the fermionic Z-f actor 



d^Zf 



=s' r — 

Jq,R idko 



GfR(q-k) GbR(q) 



(20) 



k^O 



(18) 



This expression vanishes for a local boson (Zt, = 0), which re- 
flects the fact that for a momentum-independent four-fermion 
interaction the fermionic Z-factor is renormalized only at the 
two-loop level. When inserting the particle-particle bubble 
which generates Z^ in the third line of Fig. |4l into the bosonic 
propagator in the first line of Fig. (H we observe that the flow 
equations indeed capture two-loop efl'ects. 

For the control parameter, we evaluate the diagrams with 
two external bosonic legs, obtaining the two contributions 

= g^ f GfR(k)GfR(-k)^'^ f GtRiq). (21) 

Jk,R ^ Jq,R 
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The fermionic contribution on the right-hand- side is positive 
leading to a reduction of d for decreasing A, whereas the 
bosonic contribution tends to increase 6. 

The flow of the bosonic Z-factor is obtained as the second 
frequency derivative of the fermionic particle-particle bubble: 



r 1 ^2 

dhZb = —jGfR(k -h q) GfR(-k) 

Jk,R ^ oq^ 



(22) 



^=0 



The bosonic tadpole diagram does not contribute here, as 
the 0^ -vertex u is taken as momentum- and frequency- 
independent. 

Finally, the bosonic self-interaction flows according to 



45' 
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r [GfR{-k)f[GfR{k)f 

Jk,R 



Jq,R 



{GtR{q)Y , 



(23) 



where the first terms generates u and the second, bosonic term 
tends to reduce u in the course of the flow. 

All frequency and momentum integrations in the above flow 
equations can be performed analytically. We now eliminate 
explicit A-dependences from the flow equations by employing 
the following scaling variables: 



d = 



A^Z, 



a¥z/V^V^ 



u = 



(24) 



Dependences on the Z-factors are absorbed by introducing the 
anomalous dimensions for fermionic and bosonic fields. 



r]b = - 



dlogZf 
J log A 
dlogZb 
J log A 



(25) 



The flow equations for the control parameter and the bosonic 
self-interaction are then obtained as 



dS 



J log A 

du 
J log A 



4(1 -hS) 



3/2 



(J - 3 2rjb) u-6f -\- 
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(26) 



The rescaled fermion-boson vertex g obeys the equation: 



J log A 



3-d 



(27) 



The equations (l26l |27]) have to be considered in conjunction 
with the fermion and boson anomalous dimensions: 



% = 4^ 



'7/ 



1(1 +r 



1 



(28) 



We will now investigate the fixed point of these equations and 
then solve them numerically. 



V. SOLUTION AT THE QUANTUM CRITICAL POINT 

A. Fixed point 

At the quantum critical point, the system is scale-invariant 
at low energies. The scaling regime is associated with a fixed 
point of the flow equations, which attracts the flow for small 
A. The fixed point is defined by the condition that the right 
hand sides of the flow equations for S, u, and g vanish, such 
that 



dS 



du 



dlogA dlogA dlogA 



= 0. 



(29) 



One thus has to solve three algebraic equations together with 
the equations for the anomalous exponents, Eq. (l28l) . The 
fermion-boson vertex g and the boson self-interaction u are 
relevant couplings below three dimensions, and the flow equa- 
tions for g and u have stable non-Gaussian (g 9^ 0, ii 0) so- 
lutions with finite anomalous exponents r]f and rj^. The condi- 
directly implies a relation between the anoma- 



tion 



_dg_ _ 



lous dimensions: 



3-d 



rjb 
2 



(30) 



The values of and T]t as obtained from a numerical solution 
of the fixed point equations are plotted in Fig. [5] as a function 
of dimensionality d. 




FIG. 5: (Color online) Fermion and boson anomalous exponents at 
the QCP for 1 < J < 3. 
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At quantum criticality, the fermionic properties of the sys- 
tem cannot be described in terms of conventional Fermi Hq- 
uid theory. A finite fermion anomalous dimension entails a 
fermion propagator scaling with an anomalous power law as 

Gf(AKAoj)ocA'^f-^ (31) 

for A ^ 0, corresponding to a non-analytic frequency and 
momentum dependence of the self-energy. In particular, 

i:f(0,oj)-^oj^-'^f (32) 

at the Fermi point k = 0. In two dimensions we obtain 
1-77^ = 0.68. A breakdown of Fermi liquid theory in quan- 
tum critical interacting Fermi systems is very common."^ For 
example, at a quantum critical point associated with a d-wave 
Pomeranchuk instability in two dimensions, the fermionic 
self-energy scales as cl)^/^.^"^-^^ The fermion propagator also 
develops a small anomalous dimension in QEDs^^ and in the 
Gross-Neveu model. ^^'^^ 

Concerning the collective properties of the system, a finite 
boson anomalous dimension leads to anomalous scaling of the 
order parameter propagator 

Gb(Aq, Aoj) oc A"^'-^ (33) 

for A ^ 0, where 2 - jjt, = 1.63 in two dimensions. 

In a diff'erent context, at the antiferromagnetic QCP of the 
spin-fermion model in two dimensions, Abanov et al^'^'^^ ob- 
tained an anomalous momentum scaling oc Iql"^'^^ for the spin 
susceptibility within a perturbative 1 /N calculation, where N 
is the number of hot spots on the Fermi surface. 

B. Quantum critical flows in two dimensions 

We now establish a continuous link between the micro- 
scopic bare action in Eq. ^ and the low energy behavior of 
the effective action at the QCP in two dimensions. To this 
end, we compute the renormalization group flow as a function 
of A by solving the flow equations (1261128b numerically. The 
initial conditions of our parameters are chosen to match the 
bare action in Eq. ([3]): u = 0, Zi, = 0, Zf = I, and g = I. 
The Fermi velocity and the ultraviolet cutofl' are set to unity, 
Ao = V/ = 1. To reach the quantum critical state, the initial 
value of S is tuned such that ^ ^ at the end of the flow. 

In Fig. [61(a), we show the flow of the control parameter at 
the QCP versus cutofl' scale in a double logarithmic plot. We 
observe scaling behavior 

S ~ A^-"^' (34) 

over approximately 20 (!) orders of magnitude (note the small 
values of S on the vertical axis). The scaling stops at large 
s due to tiny deviations of the tuned initial S from the criti- 
cal point. The power law scaling of S corresponds to the fixed 
point plateaus of the reseated control parameter S in Fig.[6](b), 
where also the flow of the rescaled fermion-boson vertex is 
shown. The point we make here is that, starting from the mi- 
croscopic model, the efl'ective action is really attracted by the 
quantum critical fixed point. 



(a) 




(b) 
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FIG. 6: (Color online) (a): Flow of the control parameter 6 at the 
QCP. (b): Flows of the rescaled control parameter and fermion-boson 
vertex as a function of s = -ln[A/Ao]. The ultraviolet (infrared) 
regime is on the left (right) side of the plots. 

In Fig. [7] we present flows for the fermionic and bosonic 
Z-factors Z/ and Z^ as well as the associated anomalous di- 
mensions T]f and r]h. In Fig. [7] (a), we observe that Z^ is lifted 
from zero for small s (large A) and subsequently diverges as 
a power law: Z^ ~ A"^\ which corresponds to fixed point 
plateaus of r]iy shown in Fig.[7](b). The fermionic Zf is ini- 
tially equal to unity and then also diverges as a power law, 
Zf ~ A"^^, with a slightly smaller slope than Z^, as reflected 
by the lower fixed point plateau of /// in Fig.[7](b). The nu- 
merical solution for the Z-factors fulfils the relation Eq. ([30b 
and the fixed point values match those of Fig. [5] at d = 2. 

We show flows of the -coupling in Fig.[8l Fermion fluc- 
tuations quickly generate u as is observed in the peak for small 
s in Fig. [8] (a), and then the interplay of fermionic and bosonic 
fluctuations leads to the power law scaling behavior 

accompanied by the fixed point plateaus of u in Fig.[8](b). The 
relatively fast compensation of the fermionic contribution to 
the flow of u by bosonic fluctuations, already for small s, is 
caused by the large initial bosonic propagator, G^"^° = - 1 
where S = 6c ^ 0.009 is rather small. 

Note that the rescaled variables S, g, and u start from infinite 
initial values at ^ = 0, which is simply due to the initial value 
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FIG. 7: (Color online) (a): Flows of the bosonic and fermionic 
frequency- and momentum factors Zf and Zb. (b): Scaling plateaus 
of the fermion and boson anomalous dimensions at quantum critical- 
ity versus 5 = - In [A/Aq]. r/f starts off at zero (inset) and becomes 
finite as soon as ^ 0, see text below Eq. ( [20l ). 



FIG. 8: (Color online) (a): Flow of the quartic coupling u. In- 
set: Small- 5 behavior where fermionic fluctuations start generat- 
ing u. (b): Fixed point plateau of the rescaled ^'^-coupling versus 
^= - ln[A/Ao]. 



of the bosonic Z-factor, Z^ = 0. The corresponding unsealed 
variables are finite and the flow is therefore well defined. 

Our RG equations not only yield the scaling exponents at 
the QCP, but also determine non-universal properties of the 
system. For example, the characteristic scale Aqc at which 
the quantum critical asymptotics sets in can be obtained. From 
Fig. [71 we read oflT that the anomalous dimensions attain their 
fixed point values at ^ ^ 5, and therefore Aqc ~ Aq^"^. 

Furthermore, we can compute the interaction strength Uc at 
which the system becomes quantum critical and compare it to 
the mean-field value derived in Eq. (|7]), where we have chosen 
the same values for the band cutoflT (Aq = 1) and the Fermi 
velocity (v/ = 1). Uc is the inverse of the initial value of the 
control parameter S^"^^ tuned such that S vanishes at the end 
of the flow. We find for the ratio between mean-field and RG 
interaction strengths: 



C. Quantum critical exponents 

When approaching the quantum critical point along the 
control parameter axis, the susceptibility and the correlation 
length diverge as a power law, ;^ ~ (S-Sc)~^ and^ ~ (S-Sc)~^, 
where S is the bare bosonic mass. 

The inverse susceptibility is given by the renormalized 
bosonic mass at the end of the flow,;^"^ = limA^o ^- Reading 
oflT the slope of a double-logarithmic plot of the susceptibility 
as a function of ^ - see Fig.O we find in two dimensions: 

r=1.3. (37) 

The correlation length exponent v now follows from the gen- 
eral scaling relation y = y(2 - t]i,)?'^ With rj^ = 037 in d = 2, 
one obtains 

y = 0.8 . (38) 



Uc 
IJMVT 



17.6, 



(36) 



VI. CONCLUSION 



therewith shifting the mean-field position of the QCP in Fig.[T] 
to the left. Fluctuations drastically reduce the size of the su- 
perfluid phase. 



We have analyzed the quantum phase transition between a 
semimetal and a superfluid in a model of attractively interact- 
ing fermions with a Dirac cone dispersion. The model is a 
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FIG. 9: (Color online) Logarithmic plot of the inverse susceptibility 
versus 6 - dr 'md - 1. 



prototype for systems for which the Hertz-Millis approach to 
quantum phase transitions in itinerant Fermi systems is not ap- 
plicable, since integrating out the fermions leads to a singular 
Landau-Ginzburg functional. Using the functional RG frame- 
work, we have derived coupled flow equations for the gapless 
fermionic degrees of freedom and the bosonic order parameter 
fluctuations. Fermions and bosons acquire anomalous dimen- 
sions at the QCP in dimensions J < 3. Consequently, both the 
fermion and the order parameter propagators are non-analytic 
functions of frequency and momentum. In two dimensions, 
the fermionic self-energy scales as cxp-^^ at low frequencies, 
implying the absence of fermionic quasi-particles and non- 
Fermi liquid behavior. We have also computed the suscepti- 



bility and the correlation length exponents when approaching 
the QCP along the control parameter axis. 

The strong impact of the gapless fermions on the order pa- 
rameter interaction spoils the Hertz-Millis approach to quan- 
tum phase transitions, which is based on the Landau-Ginzburg 
paradigm. The fermionic anomalous dimension 77/ influ- 
ences critical exponents in a way that resembles the inter- 
play of fermions and bosons in the Gross-Neveu model. Re- 
cently, the Gross-Neveu model has been proposed as an ef- 
fective field theory describing the quantum phase transition 
between a semimetal and an antiferromagnetic insulator for 
repulsively interacting electrons on the two-dimensional hon- 
eycomb lattice. 

In the present work we have focussed on ground state 
properties, but an extension to finite temperatures should be 
straightforward. In particular, the quantum critical regime 
above the QCP in the {6, r)-plane could be studied with the 
same truncation of the functional RG hierarchy. Another 
interesting extension would be the inclusion of the pseudo 
spin structure necessary to describe the Dirac dispersion of 
fermions on the honeycomb lattice. 

Our work may also serve as a guideline for the study of 
more complex quantum phase transitions in itinerant systems, 
for which the Hertz-Millis approach is not applicable, or at 
least questionable, such as certain magnetic transitions in low 
dimensions. 
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